Temperature of a trapped unitary Fermi gas at finite entropy 
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We present theoretical predictions for the equation of state of a harmonically-trapped Fermi gas 
in the unitary limit. Our calculations compare Monte-Carlo results with the equation of state 
of a uniform gas using three distinct perturbation schemes. We show that in experiments the 
temperature can be usefully calibrated by making use of the entropy, which is invariant during an 
adiabatic conversion into the weakly- interacting limit of molecular BEC. We predict the entropy 
dependence of the equation of state. 
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I. INTRODUCTION 

The past two years have witnessed many exciting de- 
velopments in experiments on dilute Fermi gases of ul- 
tracold 6 Li and 40 K atoms 0, S S 0, S @, 0, 01 ■ In these 
systems the interaction strength between atoms, repre- 
sented dimensionlessly as fcj?o, where kp is the Fermi 
wave-vector and a is s-wave scattering length, can be 
varied arbitrarily via a Feshbach resonance. One can 
therefore access the strongly interacting regime of — 1 < 
l/kpa < 1, where a smooth crossover from a Bardeen- 
Cooper-Schrieffer (BCS) superfluidity to a Bose-Einstein 
condensate (BEC) occurs. Of particular interest is the 
unitary limit with negligible interaction range and diver- 
gent scattering length ffc go ~ oo), in which a universal 
behavior is expected jlJ, llfj- This remarkable feature 
renders the unitary Fermi gas an intriguing many-body 
system. 

Evidence for the onset of fermionic superfluidity at the 
BCS-BEC crossover has been found in several ground- 
breaking experiments. These have measured conden- 
sate formation Q, density and momentum distribu- 
tions (3] , collective excitations Q, RF spectroscopy 0, 
and vortices Q . A thermodynamic measurement of the 
heat capacity of 6 Li atoms has also been performed very 
close to the unitary limit showing an abrupt jump 
at an estimated temperature of about 0.27Tp, where Tp 
is the Fermi temperature. However, there is no model- 
independent method to measure the temperature of a 
strongly attractive, deeply degenerate Fermi gas. In ex- 
periments the gas is characterized indirectly by an empir- 
ical temperature (or thermometry) obtained by fitting in- 
tegrated one-dimensional (ID) density profiles to an ideal 
Thomas-Fermi (TF) distribution 0. For these isolated 
systems, the entropy is a more readily measurable than 
the temperature. The goal of this work is to predict the 
equation of state in terms of measurable quantities like 
the entropy. 

In contrast to these rapid experimental advantages, 
theoretical prog ress on the crossover is quite limited 
El El El El El El El El EH- In the strongly 
correlated unitary limit, it is extremely difficult to con- 
struct a quantitative theory in terms of a well-defined 



small parameter, especially at finite temperature. There- 
fore, current theoretical understanding of experimental 
data relies heavily on the simplest BCS mean field pic- 
ture d I2I I2J. In particular, the heat capacity mea- 
surement has been explained by a crossover theory @], 
where a BCS-like ground state is generalized to accom- 
modate thermal bosonic degrees of freedom in the long 
wavelength limit f2^ |. 

In the unitary limit of a Fermi gas there are strong 
pair fluctuations in the T-matrix approximation beyond 
the mean-field level. These must be included to predict 
the energy as a function of entropy in the unitary regime. 
The entropy in turn is a function of more readily mea- 
sured temperatures in well-understood weakly interact- 
ing regimes, which are accessible via adiabatic changes 
in the coupling constant. The combination of measured 
temperatures of molecular BEC together with the state- 
equation, also allows one to estimate the corresponding 
strongly-interacting temperatures. While entropy invari- 
ance under adiabatic passage is well-understood, the cru- 
cial step here is to obtain a reliable state-equation be- 
low threshold that covers the entire range from weak to 
strong interactions. 

Without a controllable small parameter in the unitary 
regime, the determination of the equation of state is by 
no means a trivial task. In general, there are a number 
of alternative T-matrix schemes that can be implemented 
to tackle the crossover problem, and there is no a priori 
judgement of which one is the most accurate. In this re- 
gard, it is of particular relevance that Monte Carlo simu- 
lations have been performed in the unitary limit at finite 
temperature j2^|, which present useful benchmarks on 
the validity of different approximations. Here we present 
a comparative study of the equation of state of a uniform 
gas using three candidate T-matrix schemes, and show 
that, apart from a small region around the transition 
temperature, a below threshold version of the Nozieres 
and Schmitt-Rink (NSR) scheme seems to be the opti- 
mal choice for the calculations EH °f the type required 
for thermometry. We then incorporate the effect of an 
harmonic trap using a local density approximation that 
holds for large particle number. 

We characterize the entropy and temperature of the 
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strongly interacting system by an isentropic sweep to the 
weakly interacting molecular BEC regime, and a conse- 
quent determination 0, 0, 0] of the entropy S from the 
final temperature T". The T" dependence of the equa- 
tion of state is predicted. Alternatively, we can view 
this as a predicted equation of state in terms of the 
entropy, which can be readily compared to completely 
model-independent experimental observations of quanti- 
ties like the system energy as a function of entropy. 

The manuscript is organized such that the equation of 
state of a uniform unitary Fermi gas is first determined 
by comparing three different perturbation schemes, and 
then used to calculate the equation of state of a trapped 
Fermi gas within local density approximation. In turn 
the resulting unitary entropy is used to match the en- 
tropy of a weakly interacting Bose gas to obtain the de- 
sired isentropic thermometry. Finally, a summary and 
conclusions are given. 

II. EQUATION OF STATE OF HOMOGENEOUS 
UNITARY GASES 

To describe a gas of 6 Li atoms at a broad Feshbach 
resonance B = 834 G, we adopt a single-channel model 
HE m which the attractive interatomic interac- 

tion is characterized by a contact potential of strength 
U = 4irh 2 a/m. In the normal state above the superfluid 
transition temperature T c , the model can be treated per- 
turbatively with the inclusions of strong pair fluctuations 
in the many-body T-matrix approximation [22j, which 
leads to a two-particle propagator \ * w n) given by 

X (q, iv n ) = -x X! (k, iuj m ) Gb (q - k. iv n - iu) m ) , 

f l 

k..iuj m 

(1) 

and a corresponding self-energy of form, 

E (k, iu m ) = ~ J2 l + Ux{^ n f c (q " k ' Wn ~ lLUm) ■ 

(2) 

where ui m = (2m + 1) 7r//3 and v n = 2nir/f3 are, as 
usual, the fermionic and bosonic Matsubara frequen- 
cies with inverse temperature (3 — 1/fcsT. The super- 
script a, b, and c in the above two equations can be 
set to "0", indicating a non-interacting Green function 
Go(k,iw m ) = l/[iu> m — H k /2m + (j]. If this superscript 
is absent, we refer to a fully dressed (interacting) Green 
function, and thereby a Dyson equation, G (k, iuj m ) — 
G (k, iuj m ) /[l - G (k, iuj m ) £ (k, iu m )}, is required to 
self-consistently determine G and £. The only free pa- 
rameter is the chemical potential \x, which is fixed by 
using the number equation n = 2//3^ kiw G(k, iu> m )- 
According to the different combination of subscripts of 
a, b, and c, on general grounds we may expect three ob- 
vious choices of the T-matrix approximation, for which 
a nomenclature of (G Q G(,)G C will be used. These al- 
ternative choices distinguish themselves by the level of 



self-consistency contained in the Green function. 

The simplest choice, (GoGo)Go, was proposed in a 
seminal paper 13] by NSR, although these authors ap- 
proximated the Dyson equation using a leading order 
series, i.e., G — Gq + GoEGo. This fully non-self- 
consistent scheme includes the least Feynman diagrams. 
In the other extreme limit, one can consider a fully self- 
consistent T-matrix approximation of (GG) G, which is a 
so-called conserving approximation [jjj. An alternative 
intermediate scheme is to use an asymmetric form for 
the two-particle propagator, i.e., (GGo) Go H3- These 
candidate versions of T-matrix approximation have been 
extensively discussed for the attractive Hubbard model in 
the context of high-T c superconductors. The applicabil- 
ity of these theories to the single-channel model is under 
active debate, due to the complexity of numerical cal- 
culations. In this work, we have resolved the numerical 
difficulties by using a technique of adaptive step Fourier 
transformations without additional approxima- 

tion. 

The situation in the broken- symmetry state is more 
subtle. Below T c , the denominator in Eq. (|2j develops 
a pole that signals the emergence of superfluid phases. 
To remove the instability, one may extend the different 
T-matrix approaches by using a BCS mean-field ground 
state as the starting point for a perturbation theory. For 
the(GoGo) Go scheme, this strategy has been recently 
adopted and investigated in depth by Strinati et al. [la] 
and by the present authors [ijj. The advantage of the 
latter study is that a modified number equation was em- 
ployed to obtain an (approximately) correct molecular 
scattering length, which therefore makes the theory quan- 
titatively reliable over the entire crossover regime [l7j . 
The other alternative versions of (GGo) Go and (GG) G 
along this line have not been explored yet. However, on 
physical grounds, below T c one may expect only a slight 
difference among these schemes, as the ground state is 
well stabilized by the mean-field gap. 

Figure (1) presents a comparative study of these can- 
didate T-matrix approaches for a uniform unitary gas at 
finite temperature. In the superfluid state, we use the 
NSR- like formalism implemented in Ref. jjjj. The cal- 
culated chemical potentials of the gas are compared with 
a recent Monte Carlo simulation from Ref. [13 ■ Just 
above the transition temperature T c ~ 0.2Tp, all approx- 
imations appear poor, with errors of around ±10%. In 
particular, the NSR approach suffers from an unphysical 
drop with decreasing temperature, which should be due 
to the lack of self-consistency in the two-particle prop- 
agator. At high temperatures above 0.35Tp, the NSR 
result turns out to be very accurate. This observation, 
together with the excellent agreement between the NSR 
prediction and Monte Carlo data below T c , suggest that 
the fully non-self-consistent NSR scheme is generally ac- 
curate, apart from the regime just above T c . The spuri- 
ous structure about T c in this scheme might be avoided 
by a phenomenological interpolation between the results 
at T c and at a high temperature ~ OATp. The outcomes 
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Figure 1: (color online) Chemical potential of a uniform Fermi 
gas in the unitary limit as a function of temperature, in units 
of Fermi energy €f = ti 2 k F /2m — kBTp. Lines plotted are: 
our proposed results (solid lines); NSR (Go Go) Go predictions 
in Ref. (dashed line); asymmetric (GGo) Go scheme 

(dotted line); fully self-consistent (GG) G scheme (dot-dashed 
line). All these results are contrasted to the Monte Carlo data 
from Ref. 23] (open circles). 



with this idea have been denoted in the figure by solid 
lines, and referred to later as "our NSR" results. From 
this, we are able to calculate the energy and entropy. 

We emphasize that the equation of state obtained in 
this manner is certainly not ab initio, but it should be 
quantitatively valid, provided the Monte Carlo simula- 
tion data are accurate. More reliable Monte Carlo simu- 
lations with larger lattices may change these conclusions 
slightly in future, especially in the vicinity of T c , where al- 
ready there are reports of small discrepancies between the 
simulations carried out with different lattices and algo- 
rithms. The important issue here is that while there are 
discrepancies above T c which make it difficult to decide 
which approximation is better, the below T c region which 
is crucial for finite entropy thermometry methods, shows 
excellent agreement between Monte-Carlo data and our 
(GoGo)Go scheme. 

We close this section by listing in Table I the predic- 
tions for the superfluid transition temperature that were 
obtained by the various models we have discussed. The 
chemical potential, the total energy and the entropy at 
the transition temperature are also outlined. 



III. EQUATION OF STATE OF TRAPPED 
UNITARY GASES 



With the knowledge of equation of state of a uniform 
gas, we incorporate the effects of harmonic traps by us- 
ing the local density approximation, which amounts to 
determining the global chemical potential from the local 



Methods 


t c /t f 


n(T c )/e F 


E(T c )/(Ne F ) 


S{T c )/(Nk B ) 


BCS-MF 


0.50 


0.743 


1.022 


1.92 


(GoGo)Go 


0.225 


0.340 


0.290 


0.64 


(GGo)Go 


0.178 


0.508 


N.A. 


N.A. 


(GG)G 


0.150 


0.404 


N.A. 


N.A. 


our NSR 


0.225 


0.459 


0.400 


0.91 


Monte Carlo 0.23(2) 


0.45 


0.41 


0.99 



Table I: The superfluid critical temperature for a homoge- 
neous unitary Fermi gas obtained by various methods. The 
chemical potential, total energy, and entropy at the criti- 
cal temperature are also listed. We note that all the meth- 
ods, except the (GGo)Go scheme, satisfy an exact identity 
P = 2/3E that holds in the unitary limit, where P = — fi = 
— (E — TS — fiN) is the pressure of the gas. The energy and 
entropy for the (GGo)Go and (GG)G approximations are not 
available. 



equilibrium condition, 

T T F 



ft ' Mhom 



n(r) 



T F T° (r) 



(r 2 + A 2 z 2 ) , (3) 



with a density n(r) normalized to the total number of 
atoms, J n (r) dr — N. Here /ihom is our proposed NSR 
chemical potential of a uniform gas in Fig. (1), and 
lo x = uj y — uj r and lu z — Xuj r are different frequencies 
for three axis. In this case, the energy and tempera- 
ture can be taken in units of the (non-interacting) Fermi 

1/3 

energy Ep — (3Nuj x ujyLu z ) and Tp — Ep/ks, respec- 
tively. Once rc (r) are obtained, the energy and the en- 
tropy of the gas can then be calculated straightforwardly, 
using E — J n (r) Eh om dr and S = J n (r) Sh om dr, with 
our proposed equation of state of a homogeneous gas, 
E hom [n{r),T/T F (r)] and 5 hom [n(r), T/T F (r)] . 

The temperature dependence of the energy and the 
specific heat of a trapped unitary gas are given in fig- 
ures (2) and (3), and are compared to the experimental 
data obtained by Kinast et al. in the vicinity of unitary 
limit (1/kpa ~ —0.03) 0. It is worth noting that these 
experimental data are not completely raw, since the mea- 
sured empirical temperature requires a model-dependent 
theorectical recalibration [2fll |. 

Our results agree qualitatively with experimental data. 
In particular, our predictions for the specific heat show 
an apparent superfluid transition at about T c ~ 0.26Tf, 
in excellent agreement with the experimental findings 
[3oT |. The peak structure resembles the A transition for a 
BEC. The behavior of specific heat of a unitary Fermi gas 
thereby lies between that of ideal Fermi gases and that 
of BECs. The value of both theoretical and experimen- 
tal T c need further refinement, as the detailed equation 
of state of a uniform gas around the transition temper- 
ature is blurred by the semi-empirical interpolation, and 
the experimental temperatures are not known accurately. 
Another noticeable feature of our result is that it devi- 
ates greatly from the non-interacting Fermi gas behavior 
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Figure 2: (color online) Energy of a trapped unitary gas is 
plotted as a function of temperature (thick solid line), and 
compared to experimental data (open circles with error bars). 
The comparison with respect to the prediction of an ideal 
Fermi gas (dashed line) is also shown. The thin solid line is 
an empirical power law fit to the measured energies Q, i.e., 
E(T) = E [l + 97.3(T/r F ) 3 ra ] for T < T c ~ 0.27T>, and 
E(T) = E [l + 4.98(T/7>) 1 ' 43 ] above T c . Here E is the 
energy at zero temperature and T c = 0.27Tf is the experi- 
mentally extracted superfluid transition temperature Q. 
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for all the temperatures considered. This is suggestive of 
the strongly-correlated nature of unitary gases, and is in 
sharp contrast to the prediction of mean-field crossover 
theories by Chen et al. |H2^|. 



IV. ISENTROPIC THERMOMETRY FOR 
TRAPPED UNITARY GAS 

A quantitative comparison between our theorectical re- 
sults and experimental findings is prohibited, as current 
estimates of strongly-interacting Fermi gas temperatures 
are based on empirical thermometry through density pro- 
file measurements, which are strongly dependent on the 
model used, but relatively insensitive to the actual tem- 
perature A better thermometry would be obtained 
by an adiabatic conversion of the unitary gas to the deep 
BEC limit, and a consequent measurement of BEC tem- 
perature T" of a resulting Bose gas j24[- Such techniques 
have been used in recent experiments for temperature es- 
timates 0,0], but the conversion between T" and the uni- 
tary temperature T remains unknown. Below we quan- 
tify this important relation using the isentropic condition 
S'(T') = S(T). 

To be concrete, we consider a gas of 2 x 10 5 6 Li atoms 
with a potential uj r — 1555 Hz and uj z — 127 Hz, accord- 
ing to the typical setup 0. Assuming that the gas is 
swept to a field of B = 676 G 0, we find l/k F a ~ 4.9. 
Thus, the system after sweep can be well described by 
a weakly-interacting Bose gas with molecular scatter- 
ing length a m = 0.6a. The entropy S'(T'), calculated 
within the Hartree-Fock-Bogoliubov-Popov theory [32T ]. 
has been compared with the unitary entropy S(T) in Fig. 
(4a). This calibrates the unitary temperature, given in 
Fig. (4b). Accordingly, the equation of state can also 
be determined in terms of the BEC temperature T', as 
shown in Fig. (4c). It is worth noticing that this kind 
of thermometry should be very accurate, thanks to the 
demonstrated prefect adiabatic sweep of magnetic fields 
2] and a precise determination of BEC temperature 



Figure 3: (color online) Specific heat Cv = dE/dT of a 
trapped unitary gas as a function of temperature. For com- 
parison, the predictions of an ideal Fermi gas (dashed line) 
and of an ideal Bo se g as with number of particles Nb = N/2 
(dot-dashed line) |27l |. together with the experimental ex- 
tracted specific heat (line with open circles) are also pre- 
sented. It should be noted that a precise determination of 
the specific heat from experimental data of energies is not 
permissible. The experimental specific heat shown here is 
calculated from the empirical power law fit to the measured 
energies, which is outlined in Fig. 2. The resulting sharp 
jump at T c is thereby possibly an artifact of the fit. We note 
also that the abrupt discontinuity in the specific heat of an 
ideal trapped Bose gas at T c will be rounded off quickly with 
inclusions of boson-boson interactions, see, for exmaple, the 
Fig. 3 in Ref. 



V. SUMMARY 

In the present paper we have studied the equation of 
state of a trapped unitary gas with a perturbation theory 
that includes fluctuations beyond the mean-field. Our in- 
vestigations, in contrast to mean-field theory [TlL IT2I I22T ] . 
are in excellent agreement with Monte Carlo simulations 
below threshold, although there are still uncertainties of 
order 10% around threshold. We find a qualitative agree- 
ment between our theoretical predictions and the exper- 
imental findings. A fully quantitative comparison is not 
yet possible, as the experimental temperatures are not 
directly determined. We reach the conclusion that the 
empirical thermometry adopted in the experiment can be 
improved by using adiabatic sweeps at constant entropy 
to the weakly interacting molecular superfluid regime. 
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The equation of state based on this type of thermometry 
has been predicted, allowing for experimental tests. 



We emphasize that while our NSR result for the equa- 
tion of state is very accurate at both low and high temper- 
atures, its accuracy around the critical temperature T c re- 
quires a further improvement beyond T-matrix approxi- 
mation. In terms of functional-integral methods [lj], our 
NSR approach can be identified as a Gaussian expan- 
sion of the action of Cooper-pairs about the saddle point 
(mean-field solution). Therefore, a systemmatic refine- 
ment of equation of state around T c might be achieved 
by expanding the action to fourth or higher order about 
the saddle point. We will address this issue in a later 
publication. 



Our calculated equation of state can be used to pre- 
dict other observables of experimental interest, such as 
the first sound and second sound velocities of trapped 
unitary Fermi gases at finite temperature 0,0^1' which 
soon may be measured. These sound modes, particu- 
larly the second sound, may provide us more accurate 
details about the equation of state in the unitary limit, 
and therefore discriminate various approximations that 
we have marked. This extension is under investigation 
and will be reported elsewhere. 



Figure 4: (color online) (a) The temperature as a function 
of entropy in the unitary limit (solid line) and in the deep 
BEC limit of l/k F a ~ 4.9 (dashed line), with T c = 0.518TV. 
At low temperature, our calculated BEC entropysatisfies the 
scaling law for a weakly-interacting Bose gas |23|: S'(T') cx 



(T 



I-.5/2 



(b) Calibration curve of unitary temperature versus 



weakly interacting Bose gas temperature, (c) Unitary energy 
as a function of BEC temperature T'. 
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